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^ Abstract 
-t— » 

We describe a variant of the dressing method giving alternative representation of mul- 
tidimensional nonlinear PDE as a system of Integro-Differential Equations (IDEs) for 
spectral and dressing functions. In particular, it becomes single linear Partial Differen- 
tial Equation (PDE) with potentials expressed through the field of the nonlinear PDE. 
The absence of linear overdetermined system associated with nonlinear PDE creates an 
obstacle to obtain evolution of the spectral data (or dressing functions): evolution is 
■ defined by nonlinear IDE (or PDE in particular case). As an example, we consider gener- 

alization of the dressing method applicable to integrable (2+l)-dimensional iV-wave and 
Davey-Stewartson equations. Although represented algorithm does not supply an analytic 
particular solutions, this approach may have a perspective development. 

a 

£ 1 Introduction 

Completely integrable multidimensional Partial Differential Equations (PDEs) represent at- 
tractive subject of intensive study during last decades after the paper [1]. This popularity is 
due to their remarkable mathematical properties and variety of physical applications, which 
may be found in literature. Investigation approach considered in this paper is, in some sence, 
associated with so-called S-integrable PDEs [2], i.e. nonlinear PDEs which may be "linearized" 
using special technique, such as Inverse Spectral Transform (1ST) [3, 4, 5]. It is well known 
that 1ST is not the only method to study S-integrable PDEs. One may refer to Sato Theory 
[6, 7, 8, 9], Symmetry Approach [10, 11], Dressing Method [12, 13, 14, 15] . The later, in turn, 
has several formulations: Zakharov-Shabat method [12], local Riemann problem [13], nonlocal 
Riemann and <9-problem [14, 15, 16]. 

Classical S'-integrable systems are basically (1+1)- and (2+l)-dimensional. Only special 
types of multidimensional S'-integrable examples are known, such as self-dual Yang-Mills equa- 
tions [17, 18, 19, 20, 21] and the Plebanski heavenly equation [22, 23, 24]. Recently a new 
type of multidimensional partially integrable systems have been found [25], for which integra- 
tion algorithm is based on the integral operator with nontrivial kernel, which is a variant of 
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the dressing method. This recent result encourage us to search for other improvements of the 
dressing method. 

It is well known, that dressing method has been originally developed to construct nonlinear 
PDEs together with their solutions. Variant of the dressing method suggested here does not 
allow one to find analytic solutions for nonlinear PDEs. However 

1. it gives an alternative representation of largely arbitrary nonlinear PDE as nonlinear 
system of Integro-Differential Equations (IDEs). In particular case, this system becomes 
single linear PDE where potentials are expressed through the spectral function from one 
hand and through the field of original nonlinear PDE from another hand; 

2. it relates a single linear spectral evolution equation (written for some spectral function) 
with largely arbitrary nonlinear PDE. 

This is an interesting result of the paper. However, the fact that one has single linear equa- 
tion associated with given nonlinear PDE (instead of overdetermined linear system, like in 
S'-integrable case) results in system of nonlinear IDEs (or PDEs) defining evolution of the 
dressing function, which is disadvantage of our representation. Remember that dressing func- 
tions of S-integrable PDE satisfy linear PDE. As a consequence, our (largely arbitrary) PDE 
may not be derived as compatibility conditions of linear overdetermined system. 

In some sence, similar purpose (but different approach) was sought in series of papers gen- 
eralizing known (2+1)- and (l+l)-dimensional completely integrable equations. These are 
generalization of Kadomtsev-Petviashvili equation (KP) using deformation of the classical In- 
verse Spectral Transform (1ST) [28], generalizations of Korteweg-de Vries equation (KdV) and 
Nonlinear Shrodinger equation (NLS) [29], generalization of Benjamin- Omo equation (BO) [30]. 
In these papers evolution of spectral data is defined by nonlinear nonlocal equations (spectral 
data are replaced by dressing functions in our case). 

Here we start with dressing method based on the integral equation in the form [26, 27], 
where we introduce an integral operator with different type of kernel allowing us to increase 
dimensionality of PDE. As a consequence, an arbitrary function of Xi (independent variables of 
nonlinear PDE) appears in the dressing algorithm (see function $(Ai; x) in Sec. 2.1) enforcing us 
to introduce an extra constrain in the form of largely arbitrary nonlinear IDE for $(Ai; x), see 
eq.(40). Fixing function $(Ai;x), this constrain provides possibility to write single nonlinear 
PDE for single field u expressible in terms of the dressing and spectral functions. Note that 
similar extra constrain has been introduced in [25], but arbitrary function there has quite 
different origin. 

Below we concentrate on multidimensional generalizations of dressing algorithm for (2+1)- 
dimensional iV-wave equation and Davey-Stewartson equation (DS). However, generalized ver- 
sion is applicable to largely arbitrary nonlinear PDE. 

In the next section (Sec. 2) we give general algorithm deriving nonlinear iV-wave type PDE. 
We introduce an extra constrain allowing to write single nonlinear PDE for single field. Charac- 
terization of solution space for derived nonlinear PDE is given in Sec. 2.1. Sec. 3 considers similar 
generalization of dressing method for DS. Finally we represent some conclusions in Sec. 4. 
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2 Derivation of multidimensional nonlinear TV-wave equa- 
tion 

We start with usual integral equation 

$(A) = J *(A, v\ x)U(u; x)dv = *(A, i/; x) * [/(*/; x) = * * U, (1) 

where * means integration over spectral parameter appearing in both functions. There are two 
types of parameters in this equation. First, already mentioned spectral parameters denoted by 
Greek letters A, fi, v (for instance A = (Ai, . . . , AdimA)), and, second, additional parameters 
denoted by ). These additional parameters are independent variables of 

resulting nonlinear PDE. Besides, we reserve k for scalar Fourier type parameter appearing 
in integral representations of some functions. All functions are Q x Q matrices. We always 
assume dim x = M and dim A = dim fx = dim v = M + 1, where M is dimensionality of resulting 
nonlinear PDE. 

Eq.(l) is a linear equation for the spectral function U(X;x), where operator ^(A, fx;x)* is 
required to be uniquely invertible, $ is a diagonal matrix function specified below. Integration is 
over whole space of vector spectral parameter v. Function \I/(A, fx; x) is defined by the following 
formulae introducing x-dependence: 

d Xn * a p(\,fi;x)+ (hl(\)+g2 p (fi))y a p(\,fi;x) = 3> a (\;x)B2C a/3 (ii;x), (2) 

1 < n < M. 

Here C(ix; x) is a new function, which will be characterized below; h n (X) are diagonal and g n (fx) 
are arbitrary matrix functions of argument; B n are diagonal constant matrices. Short form of 
eq.(2) reads: 

L^(A, fx)V a/3 (X, fx; x) = $ a (A; x)B n a C^{fx; x), (3) 

LU X ^)(*) = + KW+9 n a ^))(*), 1 < n < M. 

Remark that derivatives d Xi ^(X, fx; x) are not separated functions of spectral parameters, unlike 
the S'-integrable case [25, 27]. Been overdetermined system of PDEs for function \1/(A, ix;x), 
eqs. (3) imply compatibility conditions, which are following: 

L^(\^)($ a (\;x)BiC a p(tx;x)) = L^(X, fx) ($ Q (A; x)B^C a p{^ x)) , n^j. (4) 

Without loss of generality, we put j — 1, and 

B 1 = I i h 1 (X)=g 1 (X) = Q, (5) 

where / is the unit matrix. Since each term in expanded form of eqs. (4) is separated function 
of parameters A and fx, these equations with j = 1 are equivalent to two sets of equations: 

d Xn ®(X;x) + h n (X)<Z>(X;x)-d xl <Z>(X;x)B n = 0, 1 < n < M, (6) 
d Xn C(fx;x) + C ln (fx;x)-B n d Xl C(fx;x) = 0, 1< n < M, (7) 

where 

C]%{fx;x) = C a e{ix;x)g n a p{fx). (8) 
3 



Eqs.(6-8) define $ and C. We refer to functions <3>, C and \I> as dressing functions, where \I> is 
expressed in terms of <£> and C due to eqs.(2). 

Thus we have specified all functions appearing in eqs. (1) and (2). Now we demonstrate 
how linear integral equation (1) is related with appropriate multidimensional nonlinear PDE 
written for fields expressible in terms of spectral function U(X; x) and dressing functions. 

System of nonlinear equations is generated by eq.(6). Derivation is very similar to derivation 
of classical integrable equations [25, 27]. First of all, we use representation for $ as ^ * U, see 
eq.(l). Then using equations (2) for derivatives *ff Xn we end up with homogeneous equations in 
the form 

V(\,H;x)*E n (n;x) = 0, (9) 
E n (X, x) = U Xn (X, x) - U X1 (A, x)B n + U(X, x)[B n , u(x)\ - Q n {X, x), K n < M, 

where function u is related with spectral function by the formula 

u(x) = C(X,x)*U(X,x) (10) 

and functions Q n satisfy the following equations: 

V(\,fi;x) *G n (fi;x) = ^ n (X,fi;x)*U(fi;x), 1 < n < M, (11) 
^{X^-x) = * aP (X,fi;x)g^(ii). (12) 

Later, function u will be field of nonlinear PDE. 

Eqs. (11, 12) along with eq.(l) will be used in Sec. 2.1 to analyse solution space of nonlinear 
system. Inverting operator in eqs. (9) one gets 

E n (X-x) :=U Xn (X-x)-U Xl (X;x)B n + U(X;x)[B n ,u(x)]-g n (X;x) =0, 1< n < M. (13) 

In the case of classical dressing method, nonlinear integrable PDE can be received for function 
u applying C(X;x)* to (13) and using eq.(7) for C Xn , n > 1. Doing the same one gets in our 
case: 

E%(x) := u Xn {x) -u Xl (x)B n + u(x)[B n ,u(x)} =H n (x), (14) 
U n (x) = [B n , u\x)] - C ln (n; x) * U(n; x) + C(/x; x) * g n {/j; x), 1< n < M, 

where function u 1 is related with spectral function by the formula similar to eq.(10): 

u\x) = C Xl (X;x)*U(X;x). (15) 

Functions u 1 and 7i n are "intermediate" functions which will be eliminated from the final 
system of nonlinear PDEs. 

System (14) has an obvious limit to classical (2 + l)-dimensional S'-integrable iV-wave equa- 
tion. In fact, if g n = for all n, (i.e. G n (X; x) = 0, H n (x) = [B r \ w 1 (x)]), then we may eliminate 
Ui using two equations (14): E% and E^, n ^ m: 

[u Xn ,B m ] - [u Xm ,B n ] + B m u xl B n - B n u Xl B m - [[u, B m ], [u, B n ]] = 0. (16) 

This is the classical (2+l)-dimensional completely integrable iV-wave equation, which has ac- 
ceptable reduction up a = u a p, where bar means complex conjugation, see, for instance, [5]. 
System (13) with Q n = becomes linear overdetermined system for eq.(16), where U(X;x) is 
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a spectral function, i.e. eq.(16) is compatibility condition for E n and E m . This is well-known 
common feature of S'-integrable models: they may be derived both algebraically through com- 
patibility condition of overdetermined linear system and using dressing method. 

However, if g n ^ for all n, then Q n ^ and system (13) may not be considered as a linear 
overdetermined system, since it has set of spectral functions, such as U(X; x) and Q n (X; x). As 
a consequence, nonlinear eqs. (14) have extra functions H n (x) and may not be received as 
compatibility condition of the system (13) through commutation of linear operators appearing 
in (13). So, similar to [25], the only way to derive system (14) from eq.(13) is the dressing 
method. 

The derived system (14) consists of (M — 1) equations and M fields, which are u and 7i n , 
1 < n < M. In other words, it is not complete. In order to write a single nonlinear PDE for 
field u we involve another important deviation from the classical approach. 

Let us split C(A; x) into two factors: 

C a p(n] x) = G^ni, x)G 2 a p(n; x), (17) 
d Xn G\^- x) - B n d Xl G\^- x) = 0, K n < M, 

d Xn G 2 (fi;x) + G ln (^x) = 0, 1 < n < M, G 1 ^; x) = G a/3 (fi; x)g^(fi), 
d Xl G 2 ^;x)=0 

where eqs.(17b-d) appear due to the eq.(7). Multiply eq.(13) by G 2 (X;x) from the left and 
integrate over A = (A 2 , . . . , Xm+i)- One gets 

E n (X i; x) := U Xn (Xv,x)-U Xl (X 1 -,x)B n + U(X 1 -,x)[B n ,u(x)}-^ n (X 1] x)=0, (18) 
1 < n < M, 

where 

L>(Ai;x) = J G 2 (X;x)U(X;x)dX, (19) 

l/ ln (Ai;a;) = - J G 2 Xn (X;x)U(X;x)d~X = J G ln (X;x)U(X;x)dX, 

g n (X 1 ;x) = J G 2 (X;x)g n (X;x)d~X, P l (X 1 ; x) = G n (X 1 ; x) - U ln (X 1 ; x), 
1 < n < M. 

We will see in the next section that off-diagonal parts of U{X\\ x) and Q n {X\\ x) have arbitrary 
dependence on x. Thus we are able to introduce one more relation among them. For instance, 
let 

M 

J2 S ap(K^i^)-^U aP (X i; x)(Bl-Bi)) = 0, (20) 
where S l a/3 are constants. Then eq.(18) gives (a ^ f3): 

M Q 



7=1 



Sip (d x ft aP {Xu x) - d xl U a p{Xu x)Bl + U ai (Xu x)u^(x)(B^ - B^)- (21) 

X 1 U a p(X 1 ;x)(B i a -B i f} )) = 0, 

M 

T, S U Bi a~ B e) = - ( 22 ) 



i=2 
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This equation is a linear equation for the spectral function U°f; additional relation (22) is 
introduced to eliminate diagonal part of U from the nonlinear term of (21). 

Multiply this equation by G^X^x) from the left, integrate over Ai and assume that 
G 1 xl (\ 1 ;x)=\ 1 G 1 (\ 1 ;x): 

M Q 

^2S i a p(dx i u a p(x)-d xl u a p(x)B i (3 + u ai {x)u lP {x){B\ - = 0, a ^ (3 (23) 

1=2 7=1 

7?W/3 

which becomes iV-wave equation if, along with (22), one requires 

Sa? = Si*i U ^ = u a ^. (24) 

Thus, nonlinear eq.(23) is equivalent to linear eq.(21) where spectral function [7°^(Ai;x) is 
related with dressing functions by the eqs. (1-8, 17,19). Detailed discussion of this relation is 
represented in the next subsection. 

2.1 Analysis of the system (1-8,17,19,21) 

In this section we characterize solution space of nonlinear equation (23) in terms of dressing 
functions $ and C. First step is solving equations (2,6,7) for \I/(A, fi; x), $(A; x) and C(/i; x). 
Eq.(2) is nonhomogeneous equation for \I/(A, //; x), so we take the following solution: 

M 

e{fi;x) J + d aj3 d{\ - n)t 

M+l 



* Q/3 (A,/i;x) = cT 1 ($ Q (A; x)C a/9 (//; x)J + <^(A - /i)e - 2 V ' , (25) 

M+l 

S(X-fi)= n^-zii) 



1=1 



(remember that dimension of spectral parameters is M + l), where is Kronecker delta 
symbol, first term is a particular solution of nonhomogeneous equation, while the second term 
is particular solution of homogeneous equation associated with eq.(2). Function (25) is not 
general solution of (2), but this is enough for our algorithm. 
Solutions of eqs. (6,7) in view of (17) read 



M 



<MA;x) = [$° a (\,k)e K *^dk, K*(\,k;x) = kx 1 + J2{kBi-hi(\))x j , (26) 

J j=2 

C a/3 (ii;x) = G 1 a (ii 1 ;x)G 2 a/3 (ii;x), (27) 
Gli^-x) = e K ^^\ Gl^;x)=e K ^C° a ^), 

/ M \ M 

Kf (/i i; x) = ii\ [x 1 + B >i , K %(V\ x) = - 9i^) x . 



31 



i=2 J j=2 



where parameter k is scalar. 
Hereafter we take 



$°(A, k) — 8(\ 2 - k)I. (28) 
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Thus expression (25) may be written in explicit form: 

K*^^% 1 ;x) + ^WqO ( ) _ f (h J a (X)+g J a g(n)\ 

<MA,ji;s) = r— 2^ + ^(A- M ) e ^V ;. (29 ) 

A 2 + /ii 

- E UiW+siaWjXj 

Due to the last term in eq.(29), eq.(l) has term e j=2 v y U a p(X]x). However, we 

would like to eliminate factor ahead of U in this term for convenience of subsequent construc- 
ts UiW+gLWjxj 

tions. To do this, we multiply eqs.(l,ll) by e J - 2 v ' : 

E u (X;x) := U(X;x) = -d" 1 (V(A; ^C^; x)) *U(fi;x) + $ 1 (\;x), (30) 
E Gn (X;x) := ^(A;x) = -a- i 1 ($ 1 (A;x)C( M ; a ;))*^( /U ;x)+ (31) 
a- 1 (V(A; x)C lB (/i; x)) * [/(/.; x) + U n (X;x), n > 1, 



where 



M 

$i(A;x) = e^= 2V y $ a (A;x) = e^ 1 (A;x) , (32) 



M 



*f(A; x) = A 2 x x + ^ {X 2 Bi + ^(A)) 
i=2 

U^(X-x) = gl a {X)U aP {X-x). (33) 



Below we need function 



G%(\;x) = GU\;x)<fo(\). (34) 



Applying J dXG 2 (X;x)- to eqs (30,31) and J dXG ln (X; x) ■ to (30) one gets equations for U, Q n 
and U n : 

U(X 1 ;x) = -Jd; i 1 (^(X 1 ]x)G\^x)^U( f i l] x)dfi 1 + 6(X 1 -,x), (35) 

^(X.-x) = -y'9- 1 (d(Ai;x)G 1 (Aii;x))(^ n (Aii;x)-^ ln (Aii;x))dAii+ (36) 

L> 2n (A i; x), 1 < n < M, 
t)" ln (Ai;x) = - J d; i 1 (^ ln (X 1 ;x)G\ii 1 -,x)^U( f x 1 -,x)dii 1 + ^ ln (X 1 -,x), 1 < n < M, (37) 



where 



$(Ai;x) = J G 2 (X;x)^ 1 (X;x)d~X, $ lri (Ai;x) = J G ln {X;x)^ 1 (X;x)d~X, (38) 
f> 2 "(A i; x) = J G 2 {X-x)U n {X-x)dX = J G 2n (X;x)U(X;x)dX, 



G 2 a }(X;x) = G 2 a ,(X;x)g^(X). 



Equation for U 2n follows from eq.(30) after applying J dXG 2n (X; x)-\ 

f> 2 "(A i; x) = - J d~ 1 1 (^ 2n (X 1 ;x)G\ii 1 ;x^U(fi 1 ;x)dfi 1 + ^ 2n (X 1 ;x), n > 1, (39) 

$ 2n (X 1 ;x) = J G 2n {X;x)^ 1 (X;x)dX. 

By construction, function $(Ai; x) has arbitrary dependence on variables x, if, for instance, 
g l a p{X) = Xi + ig l a/3 , where g l a/3 are constants, i = 2, . . . , M. Due to this fact $(Ai; x) may solve 
equation (20). Let us transform eq.(20) substituting eqs. (35-39): 

M 

E ^{^^(Ax; x) - d x X^u x)Bi - xA a p(Xi; x)(B l a - £?*)- (40) 

7=1 

9" 1 ((^^(Ai; x) - d^iXv, x)Bt ( )G 1 1 (^ x)) c/" 7/3 (/i; x)- 

where a ^ (3 and eqs. (35-39) give us 

f n (X l]X ) = - J d^(^{X l] x)G 1 ^ 1 -x))P n {ii l -x)d l i l - (41) 

J (^x„(Ai; x) - 6 Xl (Ai; x^G^in; x)U{^ x)d^ + 

^„(A 1 ;x)-6 iei (A 1 ;x)B n , K n < M. 
Deriving eqs. (40,41), we took into account an obvious relation 

- *S = " M**^, 1 < n < M. (42) 

Note, that diagonal elements of l>, 

^ a (A i; x) = | <^(A)e*( xl+E ^* x 0dA, (43) 

may be arbitrary functions of single independent variable and & 2 ™ a — <E>„™ = 0. 

System (35,40,41) represent a complete nonlinear system of equations allowing to find 
U(Xi',x) and $(Ai;x). Since u{x) = J G 1 (Xi,x)U(Xi,x)dXi, this system is alternative form 
of the nonlinear equation (23). In particular case C°(/x) = S(fXi)C°(fl), one has G^Ojx) = /, 
and this system reduces to PDE for ip(\i, x) = d~*&(Xi, x) (below a ^ (3): 

U(Xi,x) = -(p(Xi, x)u(x) + d Xl ip(\i,x), u(x) = U(0;x), (44) 

M 

E S U{ d ^ (^i^(Ai; x) - ^^(Ai; x)B i /3 - A lV 9 a/3 (A i; - flj)) - (45) 

i=2 

Q r 

E [vor7(Ai;a;)^(0;x)+ 

7=1 

9 x ^ a7 (Ai; x) - 9 Xl v? Q7 (Ai; x)-B 7 )w 7/3 (x) - Ai</? Q7 (Ai; x)u^{x){B l a - Fg)J | = 0, 



P n (X 1 -x) = -ip(\ 1 ;x)P n {Q-x)- (46) 
(vx„(Ai;a;) - ^(Ai; x)B n ^u(x) + ^^(Ai; x) - ip xlxl (X 1 ;x)B n , 
1 < n < M, 

Eqs. (44,46) with Ai = give us 

= (l + (p(0;xf) 1 (p Xl (0',x), (47) 

^ n (0;x) = (/ + ^(0;x)) '[^^(O;^-^^;^^- 

(^„(0;x) -^ 1 (0;x)B n )u(a;) 

We see that eq.(45) is linear PDE for <f°^(X 1 ; x), Ai 7^ 0, with "boundary" function (p(0;x) 

satisfying (47a). By construction, if Ai = 0, then eq.(45) is projected into (23), i.e. calculation 

of evolution of ip * (0; x) is equivalent to solving original nonlinear PDE (23). However, from 

another point of view, this evolution may be found as lim ip°* (Ai; x). 

\\ — >o 

The simplest algorithm for numerical construction of particular solutions to (23) is following. 
For given arbitrary ip(\i, x)\ XM=0 we find u(x)\ XM=0 and g n (0; x)\ XM=0 - U ln (0; x)\ XM=0 
using (47). Then solve (45) for (Ai; x)\ Xm= q. Using Tailor formulae we approximate 

tp° f (Xi;x)\ XM=A t ■ 

y? o/ (A i; x)\ XM=At » v^ o/ (Ai; x)\ XM=0 + At y^ M (A i; x)\ XM=0 . (48) 

Evolution of diagonal elements (p aa (^i',x) is fixed by y a a(Ai; x)| XM= o due to (43). Substitute 
this result into (47) we find u(x)\ XM=A t and G n (0;x)\ XM=A t - U ln (0; x)\ XM=At - Then eq.(45) 
gives Vxm(Ai; a; )|xM=At) an d so on. Solving the Initial Value Problem (IVP) (i.e. construction 
of u°f(x) for given initial data 1^(2;) ^=0) is more complicated and will not be considered 
here, since it seems to be not simpler then direct numerical solving of IVP for (23). 

Let us remark in the end of this section, that eq.(20) is not the only admittable constrain. 
Instead of zero in the rhs of this equation one might use expression L(U of {X 1 ;x);u°f{x)) which 
is linear differential operator applied to U°f(\i; x). Coefficients of this operator depend on field 
u°f(x) and its derivatives. Then expression L(U°f(\i; x); u°f(x)) appears in the rhs of (21). 
The only requirement to L is that after multiplying eq.(21) by G^(Ai;a;) from the left and 
integrating over \\ one gets nonlinear PDE for u°f . This new PDE (which replaces eq.(23)) 
may be largely arbitrary nonlinear PDE for u°* . So, as for now, represented multidimensional 
version of the dressing method is not the method for solving of nonlinear PDE, but it gives a 
new representation of nonlinear PDE. This situation is equivalent to the situation appearing 
when Fourier method is applied to PDE other then linear PDE with constant coefficients. 



3 Derivation of multidimensional Nonlinear Shrodinger 
Equation 

In the previous section we demonstrated that (largely) arbitrary nonlinear PDE can be trans- 
formed using a variant of multidimensional generalization of the dressing method for (2+1)- 
dimensional A^-wave equation. In this section we show that similar construction may be per- 
formed starting with the dressing method for (2+l)-dimensional DS. We use notations of the 
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Sec. 2. For simplicity, we take Q — 2, i.e. consider 2x2 matrix equations. Variables are 
introduced by the following system: 

d Xn * afi (\,n;x) + (hl(\)+g2 fi (fJt) > j^ a p(X,fi\x) = $ a (A; x)B2C a p(fji; x), (49) 

1 < n < M 

a XM * Q/3 (A,//;x)+(C(A)+^(//))* a/J (A,//;x) = ^$ a (A; x)Sf C a ^(/x; *) - 

^ a (X;x)B^d Xl C afi ( f i; x), 

where the first equation is identical to (2). Since Q = 2, only two B % are linearly independent, 
so we may put B % = 0, i > 2 without loss of generality. Let, in addition, B 1 = I, B M = B 2 = 
diag(l, —1), h 1 = g 1 = 0. Compatibility of (49) results in (compare with Sec. (2)): 

d X2 $(\; x) + /i 2 (A)$(A; x) - ^(A; x)B 2 = 0, (50) 
d Xn <S>(\;x) + h n (X)$(X;x) = 0, 2<n<M 
d XM $(\; x) + h M (X)<5>(X; x) - d 2 xl $(X; x)B 2 = 0, 

d X2 C(fi; x) + C 12 (^; x) - B 2 d Xl C(fi; x) = 0, (51) 
d Xn C(n\ x) + C ln (/j; x)=0, 2 < n < M, 



where 



d XM C(fi; x) + C 1M (v; x) + B 2 d 2 x C{ m x) = 0, 



C$(fi; x) = C a p(fi; x)g^(fi), I <n< M. (52) 



Eqs. (50-52) define $ and C. 

System of nonlinear equations is generated by eq.(50). Derivation is very similar to deriva- 
tion carried out in Sec. 2. First of all, we use representation for $ as \I> * U, see eq.(l). Then 
using equations (49) for derivatives *f? Xn and inverting \I>* we end up with system of linear 
equations in the form 

E 2 (X;x) := U X2 (X;x)-U Xl (X;x)B 2 + U(X;x)[B 2 ,u(x)]-g 2 (\;x) = 0, (53) 
E n (X;x) := U Xn (X;x)-g n (X;x)=0, (54) 
2 < n < M 

E M (X,x) := U XM (X,x)-U XlXl (X,x)B 2 + U(X,x)(u(x)[B 2 ,u(x)}- (55) 
2u xl (x)B 2 + K,5 2 ]) + U Xl (X,x)[B 2 ,u] - g M (X;x) = 

where functions u and u 1 are related with spectral functions by the formula 

u(x) = C(\,x)*U(\,x), u^x) = C Xl (X,x) * U(X,x) (56) 

and functions Q n satisfy the following equations: 

*(A, //; x) * g n (/i; x) = * n ( A, //; x) *£/(//; x), (57) 
^p(X,fji;x) = ^ aP (X,fji;x)g2p{fi), 1 < n < M. 

Later, function u will be field in the nonlinear PDE. 
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In the case of classical dressing method, nonlinear integrable PDE can be received for 
function u applying C(A; x)* and C X1 (A; x)* to (53), applying C(A; x)* to (55) and using eqs.(51) 
for C Xn , n > 1. Doing the same one gets in our case: 

E% 2 (x) := ^(x) -u Xl (a;)B 2 + u(a;)[S 2 ,u(a;)] = 7i 02 (a;), (58) 
E% 2 (x) := <(x) -tty^^ + u^i)^ 2 ,^!)] =H 12 (x), (59) 

£^(x) := u XAf (x) -u xlxl (x)B 2 + u(x)(u(x)[B 2 ,u(x)} - 2u Xl B 2 ^) + (60) 
ii Xl (x)[B 2 , U (a;)]=W 0M ( a ;), 

H 02 (x) = [B 2 , u^x)} - C 12 (/i; x) * U(n\ x) + C(fi; x) * £ 2 (/i; x), 
n 12 ( x ) = [B 2 ,u 2 {x)]-C 1 x 2 i {fi;x)*U{fi;x) + C Xl (^x)*g 2 {^x), 
H 0M (x) = -2u l Xi {x)B 2 - u{x)[u\x),B 2 } + u\x)[B 2 , u(x)} - [B 2 ,u 2 (x)} - 
C 1M (fi; x) * [/(//; x) + C(n\ x) * G M (fi; x), 



where function u 2 is related with spectral function by the formula similar to eq.(10): 

u 2 {x) = C xlxl (X;x)*U{X;x). (61) 

Remark that eqs.(58) coincide with (14) where n = 2 and 7i 2 = 7i 02 . Functions u 1 , u 2 and 7i m 
are "intermediate" functions which will be eliminated from the final nonlinear PDE. 

System (58-60) has an obvious limit to classical (2 + l)-dimensional S- integrable DS. In 
fact, if g 2 = g M = 0, (i.e. G 2 (X; x) = Q M (X; x) = 0, H 02 (x) = [B 2 , u\x)\, H 12 (x) = [B 2 , u 2 (x)], 
H 0M (x) = -2u 1 xl (x)B 2 -u{x)[u 1 (x),B 2 ]+u 1 (x)[B 2 ,u{x)]-[B 2 ,u 2 (x)}), then we may eliminate 
U\ and u 2 from eq.(60) using equations (58) and (59): 

£ ■= K{ 1+1 ^}-<L-<L-8u 12 u 21 u o f-A V u o f = (62) 

¥x 2 x 2 - <Px lXl = 4(ui2«2l)xixi, V? = ( u ll + «22)xi, (63) 

where 

u=( Uu Ml2 V (64) 

which is DS if % = it, i 2 = — 1, u 21 = u\ 2 - Eqs. (53) and (55) with Q 2 = Q M = become 
linear overdetermined system for this equation where spectral function is U(X;x), i.e. eq.(62) 
is compatibility condition for E 2 and Em- 

However, if g 2 ^ and g M ^ 0, then system (53,55) may not be considered as a linear 
overdetermined system, since it has set of spectral functions, such as U(X;x) and Q n (X;x). As 
a consequence, nonlinear eqs. (58-60) have extra fields H m (x), i — 0, 1, and may not be received 
as compatibility condition of the system (53,55) through commutation of linear operators. So, 
similar to Sec. 2, the only way to derive nonlinear system (58-60) from eqs. (53,55) is the dressing 
method. 

Similar to Sec. 2.1, we can take largely arbitrary equation for U°f(X-i;x) resulting to largely 
arbitrary nonlinear PDE for field u°f . For example, we want to construct such linear equation 
for U°f(Xi, x) that after multiplying it by G l (X\\ x) and integrating over Ai one gets 

M-l 

u°J M - Au°fB 2 + u of u 12 u 21 B 2 = 0, A = ^ d\ , (65) 
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which becomes multidimensional NLS if % = it, i 2 = — 1, u 2 i = u\ 2 . Let G x [\\\x) = 
A i (j 1 (Ai; x). Appropriate linear equation is following 

U°J M (Xi; x) - Af/ o/ (A i; x)B 2 + U of (X v , x)u 12 u 21 - B 2 X 2 U of (X 1 ; x)- (66) 
2( K X 1 U:{(Xr,x) + X 1 B 2 U:{(Xr,x)+XjU°f(X 1 ;x^B 2 = 

Thus nonlinear eq.(65) is equivalent to linear eq.(66) where U is expressed in terms of the 
dressing functions by the system (1,49-52). Detailed discussion of this relation is given in the 
next subsection. 

3.1 Analysis of the system (1,49-52,66 ) 

In this section we characterize solution space of nonlinear equation (65) in terms of the dressing 
functions. First step is solving equations (49-51) for \I/(A, //; a;), $(A;x) and C(/j,;x). Eqs. (49) 
represent nonhomogeneous system for ^(A,yu;x), so, similar to Sec. 2.1, we take the following 
solution: 

* Q/3 (A, 11; x) = cT 1 ($ Q (A; x)C a p(fi; x)) + 5 a p5(X - n)e *=> V / , (67) 

Solutions of eqs. (50,51) in view of (17) read 

M 

$ a (A;x) = e K ^ x \ K*(X,x)=X 2 (x 1 + x 2 B 2 a ) + x M X 2 2 B 2 a -J2hi(X)x J , (68) 

C a p(fi;x) = G 1 a (^ 1 ;x)G 2 a ^;x), (69) 
Gl^-x) = e^W), G 2 a ^;x) = e K ^C° a ^), 

M 

Kfini, x) = + x 2 B 2 a ) - x M n\Bl, K%(fi; x) = -^gipi^Xj. 

Thus expression (67) may be written in explicit form: 

*afi(\r,x) = — !^+M(A-/i)e^ 1 ; .(70) 

A2 + AH 

Equations (30-39) have the same form with 

M 

& a p(X;x) = e K ^^\ Kf{X-x) = X 2 {x 1 + x 2 B 2 ) + x M X 2 2 B 2 + Y,9U^r (71) 

i=2 

Function satisfies equation (66) where U is related with $ by eq.(35). Note, that diagonal 
elements of $ may be arbitrary functions of single independent variable, similar to Sec. 2. In 
particular case C°(/i) = S(fii)C°(jl), eq.(35) reduces to PDE (44) so that u is defined by the 
formula (47a). We see that eq.(66) in view of (44) is linear PDE for y?°^(Ai; x) with "boundary" 
function tp{0;x) satisfying (47a). Remark made in the end of Sec. 2.1 regarding numerical 
construction of particular solutions is relevant for this section as well. 
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4 Conclusions 



We applied a variant of the dressing method to derive a special representation for a largely 
arbitrary multidimensional nonlinear PDEs nonintegrable in classical sence. Although we have 
considered only iV-wave equation and NLS, reducible from the linear eqs.(21) and (66) re- 
spectively, different linear equation for the spectral function U(Xi,x) may be used. The only 
requirement is that after multiplying this equation by G 1 (Xi,x) and integrating over Ai one 
gets nonlinear PDE for u°* . 

We introduced several modifications in the classical dressing method: 

1. Eqs.(2) (or(49)) with functions h n (X) and g n ((x) showing that derivatives ^ Xj (X, fi; x) are 
not separated functions of spectral parameters. 

2. Eq.(17) splitting C(X;x). 

3. Extra constrain (20) (or (66) together with (35)) defining structure of PDE (23) (or 
(65)). This constrain is equation for function & of (Xi,x) (see for instance eqs.(40, 45) of 
the Sec. 2.1) and has no spectral origin. 

At the present form, multidimensional dressing method doesn't give explicit solutions for non- 
linear PDEs, but represents them in different form. We expect perspective development of the 
ideas outlined in this paper. 

Author thanks Prof. P.M.Santini for useful discussions of some aspects of this paper. The 
work was supported by INTAS Young Scientists Fellowship Nr. 04-83-2983, RFBR grant 04- 
01-00508 and grant NSh 1716-2003. 
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